

GENERALIZED PARITY 
RELATIONS FOR LARGE 
SPACE STRUCTURES 
WITH UNCERTAIN 
PARAMETERS 


RELATIONS FOB LARGE SPACE STBUCTUEES 
UNCEBTAIN PABAMEIEBS K.S. Ihesis 


SPACE SYSTEMS LASaPATORY 

DEPT. OF AERONAUTICS AND ASTRONAUTICS 
MASSACHUSETTS INSTITUTE OF TECHNOLOOY 
CAMBRIDBE, MA OSnSS 


(NASA-CB-176762) GENEEALIZED PABITI 



N86-2il7 30 


(Massachusetts Inst, of Tech.) 72 p 
HC A04/WF AOl < 






GENERALIZED PARITY 
RELATIONS FOR LARGE 
SPACE STRUCTURES 
WITH UNCERTAIN 
PARAMETERS 


by 

Jean R. Dutilloy 


Diplome d’Ingenieur 

de I’Ecole Nationale Superieure d’Arts et Metiers (1983) 

SUBMITTED TO THE DEPARTMENT OF 
AERONAUTICS AND ASTRONAUTICS IN PARTIAL 
FULFILLMENT OF THE REQUIREMENTS FOR 
THE DEGREE OF 

MASTER OF SCIENCE 

IN AERONAUTICS AND ASTRONAUTICS 

at the 

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 

January 1986 

Copyright © 1986 Massachusetts Institute of Technology 


Signature of Author ^ 

Department of Aeronautics and Astronautics 

January 17, 1986 


Certified by 


Wallace E. Vander Velde 
Thesis Supervisor 


Accepted by 


Eugene E. Covert 

Chairman, Department of Aeronautics and Astronautics 


GENERALIZED PARITY RELATIONS FOR LARGE SPACE 
STRUCTURES WITH UNCERTAIN PARAMETERS 

Jean R. Dutilloy 


January 1986 


2-86 


Under NASA Grant NAGl-126 


- 2 - 


Generalized Parity Relations for Large Space Structures with Uncertain 

Parameters 

by 

Jean R. Dutilloy 

Submitted to the Department of Aeronautics and 
Astronautics on January 17, 1986 in partial fulfillment of the 
requirements for the degree of Master of Science in 
Aeronautics and Astronautics. 


Abstract 


The method of generalized parity relations is one of the techniques that can be used 
to detect sensor and actuator failures on a large space structure. In this thesis, a 
model of a grid' structure is used to evaluate the performance of these relations. It 
shows their relative sensitivity to modeling errors. 

As no accurate model will be available before the structure is built in space, a 
method using sensor outputs and actuator inputs is required for the design of these 
relations. Three different estimators are studied. The second is the most 
interesting when computer memory is limited, while the third is the most accurate. 

With a few modifications, the last estimator can also generate relations optimized 
for the detection of a particular failure. They are especially interesting when the 
level of sensor noise is high. 
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Chapter 1 
Introduction 

1.1 Space Structures and Generalized Parity Relations 

Future space activities include many missions which will require the use of 
large space structures. Such structures will be lightweight and have small 
structural damping. This will allow small disturbances, in time, to produce large 
structural vibrations. For antennas and solar arrays, shape is critical, and will 
have to be controlled. Structural stress itself might be such that active shape 
control of the whole structure b required. 

Shape control will require the use of many sensors and actuators. A large 
space structure will have hundreds of these components. As they will not be 100% 
reliable, many failures can be expected during the lifetime of the structure or even 
in between periodic repairs. For example, with 200 components, each having a 
mean time between failures of 100 000 hours, we expect 17 failures in a year [3]. 
To still be able to control the vibrations, redundancy, that is more components 
than are minimally required, will be introduced. 

It is clear that the controller would have better performance if it knew which 
components have failed and stopped using them. A redundancy management 
system is designed to provide this information. It is generally made of two parts. 
First the failure detection system decides which component has failed. Then the 
reconfiguration system designs a new closed-loop configuration based on the 
remaining components. Figure 1-1 is a block diagram of the controlled structure 
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command 



Figure 1-1: Redundancy management system of a controlled space structure 
with its redundancy management system. 

Generally, failures are made detectable with the aid of a model representing 
structural vibrations. The method of generalized parity relations is one of such 
techniques. In [3] these relations were used to detect component failures on a free- 
free beam. They proved to be relatively insensitive to disturbances and sensor 
noise, but very sensitive to modeling errors. 

As these space structures are designed for a zero gee environment they need 
not support their own weight. Thus they cannot be built and tested on earth. The 
model of the vibrations will have to be based only on theoretical knowledge. This 
can give errors as high as 20% and the resulting parity relations will certainly be 
useless. This is why a technique is required that can design relations once the 
structure is in space. 
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1.2 Thesis Goal 

In this thesis, parity relations are designed to detect failures of components 
on a grid of aluminum beams. Chapter 2 presents the modeling of the grid and the 
performance of parity relations when sensor noise or modeling errors are 
introduced. Chapter 3 presents three different techniques that can build a parity 
relation using only actuator inputs and sensor outputs. In chapter 4, based on the 
same data we try to build relations with a better detection effectiveness for a 
particular failure than parity relations. 
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Chapter 2 

Parity Relations and the Space Structure 


2.1 The Parity Relation Method 


A parity relation is some linear combination of the present and passed inputs 
and outputs of a system which should be small (ideally zero) when the system is 
operating normally. 

Consider a system of L sensors and M actuators. Its linear discretized model 
based on an N dimensional state vector is 

rX(»-Hl) = ^X(0 + rU(0 (2.1) 

Y(.) = CX(0 

where i and »-l-l represent consecutive sampling times. 

^ is a N X N matrix 
ris a NxM matrix 
C is a L X N matrix 


The relations between the inputs and outputs over consecutive time steps and the 
state at time i are 


Y(0 = CX(0 

Y(i-i-i) = cx(»+i) = c^x(o + cru(o 

...and so on. 


Combining these relations over S time steps into a matrix equation, we get 
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Y(0 


C 


0 

0 

0 


u(0 

Y(f+1) 


C4> 


cr 

0 

0 


U(.-+l) 

Y(j+2) 

= 


x(0 + 

c^r 

cr 

0 


U(|--f2) 

Y(i+S-1) 


C^S -1 



... 

cr 


U(j-bS-2) 


which will be noted 

Y*(0 = C*X(0 + D*U*(0 (2.2) 

where C* is a (L S) X N matrix 

D* is a (L S)X(M (S-1)) matrix 


To form a parity relation we want to eliminate the unknown state X . This is done 
by multiplying the above relation by a vector W satisfying C* = 0 . In other 
words, is an element of the left null space of C*. Such a vector will exist only if 
S is greater than the rank of C . 

The resulting equation is 


r (0 = wV(0 - wVu*(0 = wVX(0 = 0 


(2.3) 


r(i) is called the parity residual 
is the parity vector 

will be called the parity information vector at time i 


W 


-U(0 


Equation (2.3) will hold as long as the system operates according to the linear 
model (2.1). If a sensor or an actuator fails, the relations between inputs and 
outputs are modified. Consequently the residuals based on the failed component 
will be non zero. With an appropriate set of residuals the identity of the failed 
component can be found. Of course any other perturbation such as unmodeled 
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modes, uncertan parameters or noise will also produce a non zero value. The 
effectiveness of the relation will depend on the magnitudes of these different values. 

2.2 Single Sensor and Single Actuator Relations 

The possible sets of parity relations that can be used to detect and isolate 
failures have been studied in [3]. One of particular interest is the set made of single 
sensor and single actuator parity relations. A single sensor parity relation is based 
on all actuator inputs but only on one sensor output. Similarily a single actuator 
parity relation will depend on all sensors but on only one actuator. A sensor failure 
will then have an effect on all the single actuator residuals but on only one single 
sensor residual. And an actuator failure will affect all sensor residuals but only one 
actuator residual. 

Let us consider here a system of 6 sensors and 6 actuators. Let R1 through 
R6 denote the single sensor residuals and R7 to Rl2 the single actuator residuals. 
The effect of a sensor or actuator failure on the residuals is represented in figure 
2-0, where 0 denotes a zero value and 1 a non zero value. 

Such relations can be generated as suggested in [3] by performing linear 
combinations on all the possible parity relations. If we assume C* has full rank N, 
then taking S=N-I-1 in equation (2.2) will give us a L(N+1)-N dimensional null 
space of C*. The W vectors spanning this null space generate a corresponding 
number of independent parity relations, each made of L(N-M) sensor coefficients 
and MN actuator coefficients. By combination of these relations we can eliminate 
the (L-1)(N-I-1) sensor coefficients for the single sensor relations and if M is equal or 
less than N the (M-l)N actuator coefficients for the single actuator relations. 

The main drawback of this technique is the amount of computation required. 
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Figure 2-i: Value of the residuals when a failure is present 


The number of operations needed to generate one set of relations is of the order of 
L“*N^. For the small space structure that will be studied later in this thesis, where 
we have 6 sensors , 6 actuators and 20 states, this gives about 10 million 
operations. 

In fact, single sensor parity relations can be generated very easily. Starting 
with equation (2.2) we can select the N-Pl lines involving the sensor 1 output to get 
the following : 


yi(0 


Cl 


0 0 

0 


u(0 

yi(>'+i) 




CjT 0 

0 


U(,-+l) 

yi(>'+2) 



x (0 + 

Cj4>r c^r ... 

0 


U(i+2) 

y,(i+S-l) 




... 

c]r 


U(»-l-S-2) 


where Cj is the first row of matrix C 
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which will be noted 


y/( 0 = c/x(0 -h D*u\i) 

If we assume Cj* has full rank, then taking S=N-f-l will give us one null space . 
vector Wj. The single sensor 1 relation will then be 

r,(0 = W7 Y,”(0 - Wj" Dj* U'(0 

The same can be done for the other sensors. 

It would be interesting to be able to use the same technique to generate single 
actuator relations. In fact this can be done after performing a few operations on 
equation (2.2). 

Let Dj be D without the first row of null matrices. The relation involving 
Dy*.is then 


Y(.-+l) ■ 




cr 

0 

0 


u (0 

Y(i+2) 




X(0 + 

C4>r 

cr 

0 


u(«+i) 

Y(»-l-S-l) 


C^s-I 

_ 


c^s-2r 

, 

cr 


U(j>S-2) 


which will be noted 

= c; x(o + d; u*(o (2.4) 


If we have L=M, that is if we consider only the outputs of as many sensors as 
actuators, then is a square L(S-l)xL(S-l) matrix. Provided that det(Dj,*) ^ 0, 
exbts. Multiplying relation (2.4) by gives 


u,‘(0 = -D/‘ C*, x(0 + D,-' Y/(0 


(2.5) 


Equation (2.5) is the equivalent of (2.2) with taking the place of Y and Y^ of 
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U*. We can then select the S-1 lines involving actuator 1 input, and by taking 
S=N+2 generate a single actuator 1 parity relation. 

The problem then reduces to the computation of and C^*. This 
would require an important amount of computation if D^.* were just any ordinary 
matrix. However its pseudo triangular form will be helpful here. 

First lets find when D^* can be inverted. Using the relation 

ab] 

det CD = det(A).det(D - C A'^ B) 

gives det(Dj ) = ( det(C f) )^'^. Then exists if det(C 1)^0 . 

Now find an expression for D^.'^ Assuming it is pseudo triangular like D^* we 

have 

’Xq 0 ... 0 

X, Xq ... 0 

= 

^-2 ^-3 •" ^ 

Solving Dp* = I^s-i) us the relations 

Xo(cr) = ij^ 

Xi(cr)H-Xo(c^r) = o 
X2(cr) + Xi(c^r)-i-Xo(c^2^) = o 
... and so on. 



Which leads to the solution 
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Xo = (crr‘ 

Xj = -XQC^r(cr)-i 

X 2 = - ( Xj c r+ X q c 4^2 j.) ( c rr^ 

x^ = -{x^c^ r + x^c ^^ t + XqC^^ r){ cry^ 

... and so on 


For implementation on a computer, the generation of the matrices X^ can be 
programmed recursively. Starting with 

Xo = (cr)-^ , xXq = ( c r)-^ c ^ and G=r{cry^ 


we get from step n-1 to step n with the relations 


n 

XX. 


= - ( XX^J ) G 
= ( X„ C H- XX 


n-l 




This also gives us the result of as we have 

XqC^P 


D 


c*. = 


x^c^ + x^c^ 




X^C^ + X^C^^ + XqC^ 


= 

XXo' 

XXj 

XXo 

J 

"■ 


The number of operations required by this technique is approximately 
4M^N^+2MN^. Consequently this represents a dramatic improvement for a large 
space structure where M and N are large. In the case of our small structure the 
number of operations required to generate a set of actuator relations now reduces 
to about 150000 from 10 million. 


Another advantage is that now the single sensor and actuator relations are 
generated from an identical matrix relation. In each case the W vector is made of 
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the coefficients of the single element. As we will see this is very useful when we 
want a relation to be more sensitive to a particular failure. 


2.3 Modeling of the Space Structure 

A model representing a space structure was created to study the effectiveness 
of parity relations. This model represents an experimental apparatus used at 
NASA Langley Research Center to demonstrate control techniques for large space 
structures. The structure is a grid of aluminium beams controlled by 6 inertia 
wheel actuators. Position sensors and rate gyros measure the vibrations of the 
structure in the direction perpendicular to the grid. 

The model is a state space representation based on the 10 lowest frequency 
modes and coresponding mode shapes and damping coefficients. 


Freq. w in rd/s 

2.33 

4.05 

9.04 

13.9 

19.5 

30.0 

36.0 


46.7 

65.1 

Damping coef. ^ 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 


Table 2-1: Frequencies and damping coefficients of the model 


The state vector is made of the 10 modal amplitudes and 10 modal velocities. 
01 02 ^2 - ^10 ^10 1 

The input vector b made of the torques of the 6 actuators. 

U'' = [ Uj Ug Ug U^ Ug Ug ] 

The model takes into account the output of 3 position sensors and 3 rate gyros. 

Y''= [yj y 2 Y3 74 75 76 1 
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The matrices relating X,U and Y in the continuous case will be noted 


f X(0 = AX(0 -I- BU(0 

1 Y(0 = C X{t) 


Let Wj be the frequency of mode i and its damping coefficient. Let py, gjj^ and ajj 
be the modal coefficients of mode i at position sensor location j, gyro location k and 
actuator location 1. Then as shown in [3] we have 


And 


0 

1 

0 

0 

- u \ 


0 

0 

0 

0 

0 

1 

0 

0 

" 2 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

^ 1.1 

^ 1.2 


^ 1.6 

0 

0 


0 

^ 2.1 

^ 2.2 


^ 2.6 

0 

0 


0 

^ 10.1 

• 

^ 10.2 


^ 10.6 


* 

Pl.l 

0 

P2.1 

0 

Pl .2 

0 

P2.2 

0 

Pl .3 

0 

P2.3 

0 

0 

Sl.l 

0 

S2.1 

0 

Sl.2 

0 

S2.2 

0 

Sl .3 

0 

S2.3 


0 0 
0 0 
0 0 
0 0 


0 1 
-2?,o«,o 


Pio.l 

0 

P10.2 

0 

P10.3 

0 

0 

SlO.l 

0 

S10.2 

0 

S10.3 


The use of parity relations requires that a discrete model of the structure be 
obtained. Let T be the sampling time, then our model becomes 
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x{i+i) = -h ru(o 

Y(.} = CX(0 

with ^ = e^'^ = I-hAT + A2 T^/2 -t- ... 

and r= [ Jq <f(t) dt ] B = I T -!- A T^/2 -1- A^ T^S! + ... 

2.4 Parameters for the Simulations 

The results of simulations involving parity relations will be highly dependent 
on the level of input and output signals and sensor noise. The ranges allowed on 
the experimental apparatus are 20 ounce-feet for the actuators and 1 inch for the 
position sensors. To simulate the constant movement of the structure we want to 
apply a random input to our model. But has the structure as very small damping, 
the level of torques we can apply without exceeding the range of 1 inch is very 
small. To be able to apply torques of higher level a state feedback controller- was 
introduced. The input signal was then the response of the controlled system plus a> 
small random signal. 

Sensor noise can be expected to produce an error of 10"* to 10"^ feet on the 
position sensors and 10'“* to 10'^ rd/s on the rate gyros. The two extremes were 
simulated by adding uniformly distributed errors between 10"* feet and 10"* rd/s or 
10'^ feet and 10"^ rd/s. 

Finally, we want to have a constant level of output in the simulations. To do 
this, state vectors were generated over 1000 time steps starting from X=0 . The 
last vector was then stored and used as the initial state of the simulation. 
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2.5 Simulation of Parity Relations and the Effect of Noise 

Single sensor and single actuator parity relations were generated for the 
model of the grid using the techniques described previously. They are composed 
respectively of 146 and 147 coefficients. The sampling time was taken to be 0.10 
sec. 

But we still have one degree of freedom in the choice of our relations which 
corresponds to the magnitude of the null space vector W. In order to get 
comparable relations we will require that the norm of W be such that the effect of 
sensor noise on the residuals is the same for all relations. This effect is measured 
through the covariance of the residuals when ho failure is present. Let nj be the 
noise signal of sensor i and N the noise vector coresponding to Y*. 

nIO =.( ni(0 •" ••• 1 

Our relation (2.3) becomes with noise 

t{{) = + N{i)) - WVU*(0 = W^N(0 (2.6) 

Then the covariance of the residuals is 

E(r2(0) = N(0 N^(0 )^ = Zj 

Here Wj and nj(0 are the components of W and N(i). 

Let ||W||j^ = [ E'(nj^(j)) be the N norm of W. Then choosing ||W||j^= l 

or any other constant will give us the same level of residuals in the no-fail case. 

As in our case all the £'(nj^(t)) have the same value q, || ||jyj is the euclidian 
norm. Then choosing ||W|| = 1 for all the relations will give 


(2.7) 
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E(At)) = ||W||2 q = , 

Simulations involving these relations with the model of the structure were 
carried out over 150 time steps. Two failures were introduced : 

A "zero” failure of sensor 1 (y^(i}=0) from time step 40 to 80 
A ”zero” failure of actuator 1 (uj(0=O) from time step 110 to 150 

The first simulation is the ideal case where the sensor outputs are not 
affected by noise. Thus the residuals are zero when no failure is present or if the 
failure does not affect the residual, as is indicated in figure 2-0. Four of the 
residuals are presented in figure 2-2 to 2-5. Rl and R7 are the two residuals 
affected by both failures. R2 is one of the single sensor residuals that should not be 
affected by sensor 1 failure and R8 is its equivalent for actuator 1 failure. The two 
important failure' signatures are the effect of sensor 1 failure on Rl and actuator 1 
failure on R7. These are the ones that enable us to isolate the-failed component. 
We can already see that sensor failures are easier to detect as their signature is 
much larger than the actuator one. The residuals tend to grow at the end of the 
simulation because the system response level grows after failure of actuator 1 due 
to poorer performance of the control system. 

In the second simulation, sensor noise of 10"* feet and 10"* rd/s was 
introduced. The effect on residual Rl and R7 is shown in figures 2-6 and 2-7. The 
perturbation due to this level of noise is si^l. 

Finally sensor noise of 10*^ feet and 10*^ rd/s was introduced. This 

corresponds to figure 2-8 and 2-9. Rl still allows us to detect the sensor failure. 
However for R7 the noise on the residuals covers the actuator failure signature and 
detection is impossible. 


RESIDUAL R2 

.040 -. 020 .000 .020 . 040 



Sensor 1 failure Actuator 1 failure 


Figure 2-2:, Residual R1 without sensor noise 



Sensor 1 failure Actuator 1 failure 


Figure 2-3; Residual R2 without sensor noise 
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Figure 2-5: Residual R8 without sensor noise 
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Figure 2-7: Residual R7 with sensor noise of 10'^ ft and 10“* rd/s 
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2.6 Simulation with Parameter Mismatch 

Having seen the performance of parity relations generated using the accurate 
parameters of the structure, it is now interesting to see how uncertainties in the 
knowledge of these parameters can affect the residuals. 

Uncertainties were created using one model to generate the relations and 
another slighty different for the simulation. The mismatches were created using 
the law 

A 

Pm = P ( ^ + 100 r ) 

where p^^ is the mismatched parameter 
p is the original parameter 
A is the percentage of mismatch 

r is a random number uniformely distributed between -1 and +1 

This law was applied to the frequencies and mode shape coefficients in the 
continuous case. The model was then discretized. 

The easiest way to carry out these simulations would have been to use the 
mismatched models with the original set of parity relations. However the result 
would biased by the fact that the mismatch also affects the controller. To avoid 
this a set of relations was generated for each mismatched model and used with the 
true model of the structure. 

Four models were generated. Two corresponding to a 5% mismatch and two 
for 10% . The worst case for 5% is presented in figures 2-10 to 2-13. We can see 
that with Rl the sensor 1 failure can be detected. The actuator failure however 
cannot be detected as its effect on R7 is covered by the effect of the mismatch. As 
for the worst 10% case presented in figures 2-14 to 2-17, no failure can be detected. 
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Figure 2-10; Residual Rl with 5% parameter mismatch 
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Figure 2-11: Residual R2 with 5% parameter mismatch 
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Figure 2-13: Residual R8 with 5% parameter mismatch 
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• Figure 2-16; Residual R7 with 10% parameter mismatch 



Figure 2-17; Residual R8 with 10% parameter mismatch 
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Chapter 3 

Estimation of Parity Relations 

3.1 Discussion of the Estimation Algorithm 

We have seen that an error of 5 % in the knowledge of the parameters of the 
structure will give us useless parity relations for the detection and isolation of 
actuator failures. Much greater uncertainties can be expected in the a priori 
knowledge of the parameters. Consequently a reestimation of the coefficients of 
each relation will have to be performed once the structure is built in space. 

This could be done by running an on line identification of the parameters of 
the structure and using this data to generate parity relations. Another solution is 
to directly reestimate the coefficients of the relations through an on line 
identification. Assuming that the structure will not suffer any failure at the 
beginning of its life, we will then have available the inputs and outputs in the no- 
fail case. Based on this data the estimation algorithm will try to find the relation 
that gives us the smallest residual. 

As R7, and more generally, the single actuator residuals are theimost affected 
by mismatch and noise, only the estimation of the single actuator 1 parity relation 
will be considered in this chapter. The algorithms are however valid for any parity 
relation. 


-32- 


3.2 Kalman Filter Estimation 


A parity relation can be estimated using the Kalman filter algorithm. Let the 
coefficients of the desired parity relation be the state of the noise free stationary 
process 


P,(i+1) = P,(0 (3.1) 

i being the time step of the estimation. 


Under the no-failure hypothesis this relation should give us a zero residual in the 
noise-free case, and a noisy residual if some sensor noise is considered. This 
residual can be considered as a measurement of the process (3.1). Let n be the 
noise on the residual, n = W N according to (2.6). A measurement equation can be 
written as 

y(0 = P]'Pi(0 - n = r^(0 - n 

where y is equal to 0 and n has covariance q as shown in (2.7). 


Let P^ be our estimate of the desired parity relation P^.. The Kalman filter 
equations are 





= p(0 Pi(0 1 P[(0 P-(0 Pi(0 + q ]-‘ 


Measurement ^ 
incorporation 


p+(0 


= P,-(0 + K(0(o-P»Pi(0) 

= p-(i) - K(0 p](0 p-(0 


Time 

update 


Jp;{i+i) = p/(.) 
= p+(0 


where P is the estimation error covariance and K the Kalman filter gain. 


-33- 


However these equations could lead to the trivial solution = 0. To avoid 
this we will require that the norm of the sensor coefficients of be kept equal to 1 
as in (2.7). This will be done by changing our process equation to a norm taking 
equation 




which gives the corresponding time update equations 


Time Pg*(«+1) = Pg"*"!*) / norm(i) 

update < 

p-{i+l) = P+(0 / norm2(0 

where norm(i) is the norm of the sensor coefficients of Pg(0- 


An estimation using this algorithm was carried out over 200 time steps. The 
initial parity coefficients are those of the 5% missmatch case presented in figure 
2-12. No sensor noise was introduced. However q was set to .33 10'* which 
corresponds to sensor noise of the 10"* case. Figure 3-1 shows the estimated 
residual r^=Pg(0 Pj(0 over the estimation sequence. Figure 3-2 shows the residuals 
given by the initial coefficients and the same information vector. The comparison 
of both allows us to visualize the effect of the filter. In fact it shows that the filter 
reduces the residual due to model mismatch, but makes no more progress after 50 
time steps. 

The last estimated relation was then used in a simulation identical to the 
ones performed in chapter 2. The result is given in figure 3-3. As far as the 
Kalman filter is concerned the results are what we expected. The residuals in the 

Q 

no-fail case were reduced to a level slightly higher than the .33 10* level given by 
the parameter q. The problem is that the resulting parity relation is useless as the 
actuator failure signature has almost disappeared. To get a relation that will still 
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show th«^ effect of the actuator failure we » must estimate the coefficients more 
accurately. But this can only be done if the sensor noise level is smaller than that 
considered here. 



Figure 3-3; Relation estimated with the Kalman filter 

Two other drawbacks of this technique are the computation time and amount 
of memory required by this algorithm. The number of operations per time step is 
3R^ where R is the number of coefficients of the parity relation. In our case with 
147 coefHcient this gives around 65000 operations. This is very likely to be above 
the capacity of a space based computer. Hopefully we need not use consecutive 
time steps in our estimation. We could for example take only one parity 
information vector every 10 time step thus reducing by 10 the computation power 
required. The algorithm also requires the storage of ^/ 2 N^ numbers for the 
covariance matrix P. This might be a major problem as parity relations for large 
space structures will be made of hundreds of coefficients. 
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3.3 Estimation through Minimization of a Structural Distance 

We have seen that the Kalman filter algorithm requires an important amount 
of memory which might not be available on a space based computer. This is why 
developing an algorithm requiring less storage capacity would be of some interest. 
The minimization of a structural distance is one such algorithm. This technique is 
described in [4]. The basic points will be presented here. 

Let Pg be our estimated parity relation and r^ the resulting residual. Let P^. 
and r^ be the desired relation and residual. We have 

rJO = pI( 0 Pi(0 
r,(0 = P>1 Pi(0 

where r^(i) = 0 in the no-fail case. 

We want to minimize the distance D(t) 

D(0 = [Pe(0-P,(0rP[Pe(0-P,(0] 

where P is a symmetric weighing matrix. 

Let A(j) = Pg(») - Pj(0- We have 

D(t+1) - D(0 = [A(«+l) - A(0fP (A(i+1) - A(0] + 2 A^(«) P [A(i-Hl) - A(0] 

As A(t) is unknown, we impose that it follows 

A(t-bl)-A(0 = hP-lPj(0 
where h is a scalar. 


Then we have 
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D((+l) - D(0 = h ( h pjo p-‘ p.(i) 2 Pjo Pi(0 1 

Let 

X P>-) Pi(0 

h =: 

pI(0 p-' Pj(0 

our distance will be decreasing for 0 < \ < 2 , the optimum being X=1 . The 
identification algorithm will be 

r,(0 P-‘ Pi(0 

P,(‘4-l) = P,(0 - — ^ ; 

Pj(0 P-I P,(0 

X is used to reduce the effect of noise on the estimation. For 0 < X < 1, reducing 
the value of X will reduce the effect of noise but will also increase the time required 
for the estimation. The P matrix is used to modify the importance of the different 
coefficients of in the distance. The more importance the coefficient has, the 
more quickly it will be identified. 

As in the case of the Kalman filter this estimation process would ultimately 
give us the trivial solution P^ = 0. To avoid this we will restore the norm of the 
sensor coefficients to 1 after each identification step. 

Pe(*+1) = Pe(‘+1) / norm(i+l) 

The algorithm was first tried with X=1 and P the identity matrix. In this 
case a geometric interpretation of the algorithm can be given. Consider a parity 
relation made of 3 coefficients (0,0,1), and the corresponding information vectors 
[a,b,0| where a and b can take any value. Lets assume that the first 2 coefficients 
correspond to the X and Y directions of figure 3-4 and the third to the Z direction 
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orthogonal to the figure. If our initial estimate of the relation is Pg(0)=(l,l,l] and 
we have the information vectors P[(1)=|0,1,0) and P](2)=[l,l,0| then our two 
successive estimates of P^ will be Pj(l)=[l,0,l] and Pj(2)=[l/2,-l/2,lj. The x,y 
component of these vectors are represented on figure 3-4. As shown by the figure, 
the algorithm is equivalent to an orthogonal projection of PJi) in the direction 



Figure 3-4: Successive estimates given by the structural distance 

The estimation with X=1 and P=I was carried out over 3050 time steps 
without introducing any sensor noise. The initial relation is the 5% case of figure 
2-12. Figure 3-5 shows the estimated residual over periods of 50 time steps 
covering the estimation time. Figure 3-6 shows the residuals given by the initial 
relation. The algorithm is not too succesful in reducing the covariance of the 
residuals. This can be blamed on the bad weighting of the gains. Actuator 
coefficients are much larger than their sensor counterparts and the algorithm puts 
all the weight on them. Consequently the sensor coefficients are almost kept 
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constant and the possible reduction of covariance is limited. 

We have to introduce some weighting P in the structural distance. One idea 
is to find P such that the distance is equivalent to the covariance of the residuals of 
the estimated relation P^(0- Lets suppose that P = £'(Pj P^^). Then 

D(0 = [Pe(0-Prr^(PiP[)[Pe(0-Prl • 

As Pg and P^ are not random 

D(o = -El (p,{.) - p,r(Pj p ] ) (p.(.) - p,) ) 

= El r.2(0 - 2 r,(0r, - | 

But as r=0 by definition of P, 

D(0 = El r,2(0 I 

The problem is that such a weighing would require the storage of ^/ 2 N“ coefficients 
as in the Kalman filter case. To avoid this only the diagonal terms of £'(Pj PT) were 
taken to create a diagonal matrix. 


Diagonal (P) = ( y^ , , yj , yg , - , 72 . - . Vq . - - Xe » ^1 > - > 1 

S times S times S times S times 

The result of the new algorithm is shown in figure 3-7. The residuals are now 
comparable to those achieved with the Kalman filter. The algorithm is not able to 
give a better estimate even though no noise is simulated. The usual simulation 
using the last estimated relation is presented in figure 3-8. 

Using the same weighting matrix but X=0.1, sensor noise at the 10"^ level 
was introduced. The result of the estimation is shown in figure 3-9 compared with 
the 5% mismatch residual in figure 3-10. The simulation involving the estimated 
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Figure 3-8: Relation estimated without noise using the structural distance 
relation is shown in figure 3-11. This shows that the algorithm still performs 


RESIDUFIL Ft7 ESTIMRTEO R7 

.0040 -.0020 .0000 .0020 . 0040 * -.0040 -;-0020 .0000 - 0020 .0040 


-] 1 r- 1 1 1 I I 

0 25 5<5'lOOO lozs /oso'iooo 2 025 20SO'3000 3 o 25 so^O 

TIME STEP 

Figure 3-9; Residual of the structural distance with noise 
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Figure 3-10; Residual of the 5% mismatch case 
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Figure 3-11; Relation estimated with noise using the structural distance 
correctly with sensor noise as long as X is reduced sufficiently. 

3.4 Eigenvalue Decomposition of a Parity Information Matrix 

The parity vector made of the R coefficients of a parity relation and the 
corresponding parity information vectors Pj are elements of the same R dimensional 
vector space. The definition of a parity relation states that when no noise is 
considered we have 

r(0 = PTPi(0 = O 

The vector space interpretation of this definition is that the parity vector is the 
vector orthogonal to all parity information vectors. This property was used in the 
last section when we estimated the parity relation by performing successive 


f 
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orthogonal projections of our estimate in the direction of the information vector Pj. 

Now instead of performing an estimation at each step we build a matrix A 
made of different information vectors. 

A = Pi(l) Pi(2) ... Pj(n) 

According to the definition the corresponding parity relation is in the left null space 
of this matrix. Consequently, without noise, a null space algorithm based on a 
matrix containing at least R-1 information vectors gives us the corresponding 
parity relation. 

If some noise is introduced the exact null space vector no longer exists. 
Instead we have .to find the vector ”as orthogonal as possible” to A as it is called in 
(Ij. This is the vector that minimizes the covariance of the residuals E{t^). 

This vector can be found by performing an eigenvalue decomposition of A A^. 
First we note that 

T [ 1 [ Pjl) 

A A = Pi(l) ... P.(n) 

L J 

= E ( Pi(j) pTU) I 

j=l 

Then if our estimation is based on n parity information vectors, the covariance of 
the residual r^ corresponding to the relation P^ is 
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=- El(PTPi(j))(Pi(i)Pr)l 

“ j=l 

= - pT ( i: (Pi(J) pl(j)) 1 Pr 

“ j=l 
=- p;(AA^)p, 

Q 


Now if the decomposition of A A^is 


A A = 


R 


where ^ ^.nd the Vj^ are orthonormal 

and the Xj^ are ordered by magnitude 


'R 


v; 


V, 


we have 


=7E\IpXP 

" j=i 


If we impose that the norm of P^ be one, then 

IIPJI =EIpXP = ‘ 

j=l 

And consequently E{t^) is minimized if P^, is the eigenvector corresponding to 
the smallest eigenvalue Xj. This vector is ”as orthogonal as possible” to A. 

In our case imposing the norm of P^ to be one is almost the equivalent of 
having the norm of the sensor coefficients to be one as they are much larger than 
the actuator coefficients. However if this approximation is not valid, for example if 
the sensors have different noise levels, then the algorithm can be modified as 
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follows. Let II llj^ represent the N norm that takes into account the different noise 
levels. Then E(t^^) can be written 




1 R 

° J=1 




But as the eigenvectors Vj^ are an orthogonal basis of our vector space, can be 
written 


P, =ElPTVk|Vi 

i=i 

and we have 

I|P,IIn =E IPr^kPlIVkllw"— 1 

j=l 


The best relation is given by the eigenvector Vj^ corresponding to the smallest ratio 


In [1] a singular value decomposition of A is used. The parity relation is 
given by the left singular vector corresponding to the smallest singular value. This 
technique might be more accurate for small eigenvalues (2). The eigenvalue 
technique, however is more interesting when memory capacity is limited. Storing A 
requires keeping R by n numbers. As A A^ is symmetric and can be built by adding 
up the Pj(0 P](t) matrices it only requires the storage of numbers. 


Our single actuator relation is made of 147 coefficients. The estimation 
would then in theory require 146 time steps. But if the system is slow compared to 
the sampling rate these vectors might not be representative of the information 
vector population. To be sure to get an accurate estimation some redundancy was 
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iatroduced by taking 200 time steps* 

In the first run, no sensor noise was introduced. The relation obtained was 
used in the two failures simulation and the result is shown in figure 3-12. We have 
here a true parity relation as the residuals are zero in the no fail case. The 
interesting point is that the failure signatures are quite different from what we had 
in figure 2-4, even though the inputs and outputs used in both cases are the same. 
According to section 2.2 the null space vector with S=N-|-2 is of dimension 1. 
Consequently the normalized parity relation is unique. In practice we find that 
normalized relations with very different coefficients can have almost zero residuals 
in the no fail case. 



Sensor 1 failure Actuator 1 failure 

Figure 3-12; Relation estimated without noise using the eigenvalue method 

In the second run, sensor noise at the 10’‘* level was introduced in the 
estimation and simulation. The result is shown in figure 3-13. The residuals in the 
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no fail case are comparable with what we had with the two other estimators. But 
here the failure si^ature is larger so the effectiveness of the relation is improved. 
But as usual the introduction of noise resulted in a reduction of the failure 
signature. 



Figure 3-13: Relation estimated with noise using the eigenvalue method 

Finally we note that this algorithm through the eigenvectors generates an 
orthogonal basis of all possible parity relations. For each we have a measurement 
of the residual in the no-fail case given by the corresponding eigenvalue. For the 
estimation with sensor noise the 10 smallest eigenvalues are given in table 3-1. 
They are very close to one another which means that the 10 corresponding relations 
will have almost the same residuals in the no-fail case. But their respective failure 
signatures could be very different. Figure 3-14 shows the behavior of relation 
number 5 in the two failures simulation. 
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relation number 

D 

2 

3 

4 

5 

6 

D 

8 

9 

10 

eigenvalue in 10"^ 

0.97 

1.23 

1.54 

1.76 

2.21 

2.24 

2.42 

2.59 

2.66 

2.70 


Table 3-1; 10 smallest eigenvalues of the information matrix 
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Chapter 4 

Optimization of a Covariance Ratio 

4.1 The Covariance Ratio and the Failure Signature 

In section 3.4 we saw that as sensor noise is introduced many different 
relations tend to have equivalent level of residuals in the no fail case. Up to now we 
selected the one with the lowest level so as to be as close as possible to a parity 
relation. This is not always the best choice. Comparing figures 3-13 and 3-14 we 
find that relation 5 corresponding to 3-14 is more suitable for failure detection. 
Both relations have a comparable level of residuals in the no-fail case but relation 5 
has a much larger actuator 1 failure signature. Estimation algorithms could be 
greatly improved if we could measure the performance of the relations they 
generate. A good criterion for this would be to compute the ratio of the 
covariances of the residuals in the fail and the no-fail case. To select the single 
actuator 1 relation, for example, this would be 

jE'[r^(i)] with failure of actuator 1 

= 

^[*■^(01 without failure 

Any linear relation that can be generated by our estimation algorithm is of 


the form 
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r(0 = w;Y*(0 + W^U’(0 

W, is the vector of sensor coefficients 
8 

W is the vector of actuator coefficients 
a 

r is the residual 

The covariance of this residual can be computed quite accurately in the no-fail case 
as the corresponding inputs and outputs are available. On the other hand, the 
covariance when a failure is present will have to be estimated. First we must find 
how a failure affects the residuals. Lets assume that we want to detect an actuator 

4t 

1 failure. A failure appears when the input computed by the controller Uj and the 
input applied to the system are different. The residual becomes 

r^(0 = WjY*(.l-bW^Ui*(0 

This can be written 

MO = wj y”( 0 + u,;(0 -I- |u,*(0 - u„*(.-)l 

As Ujg and Y* are the no-fail inputs and outputs of the system, the residual they 
generate, r^, is equivalent to the residual in the no-fail case, r^^. Consequently, the 
signal Uj (») - (») is the one that allows us to detect the failure. It affects the 

residual through to give the failure signature W^[Uj (i") - Uj^ (i)]. 

If we assume r has the same covariance as r„, the covariance ratio becomes 
s n’ 


- Ui3*(0) (Ui*(»l - Uj/iori 



If also the failure is such that Ujg*(0 is known, then the covariance ratio can be 
estimated based only on the inputs and outputs in the no-fail case. 
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4.2 Model Based Covariance Ratios 

To see how thb new technique performs we first assume we have available an 
accurate model of the structure. We want to find the single actuator 1 relation 
with the best covariance ratio for a zero failure of that actuator. 

Going back to relation (2.5), after selecting the lines involving the actuator 1 
input we have a relation of the form 

Ui*(0 = C* X(0 + D* Y*(0 

where C* and D* represent the selected rows of C^* and D^.*^ 

This relation can be modified to take into account sensor noise given by relation 
(2.6). We get 

Uj*(0 = C* X(0 + D* [Y*(0 + N{t-)1 ■ (4.1) 

Previously, to build a parity relation, we multiplied this relation by the vector 
W satisfying W^ C* = 0. To have such a vector the number S of coefficients in W 
had to be at least one greater than the rank N of C*. Now we will choose this 
vector, renamed W^, so as to optimize the covariance ratio and get what will be 
called an optimized relation. This means that W^ can be of any size. However if the 
dimension of W^^ is greater than N, as we will show later, we are guaranteed that 
our relation will perform better than a parity relation. Multiplying by W,, relation 
(4.1) becomes 

wj u,’(0 = c* X(0 + wj d' y*(0 + wJd' n(0 


and the corresponding residual is 
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= - < D* Y*(0 + u*{t) = C* X(0 + w; D* N(t-) 

Now the unknown state will be an added source of noise. What we hope is that 
this will be largely compensated by an increase in the failure signature. 

The covariance of this residual is given by 

«If„^(0l = w^c*£TX(0x"(.)l c*X + <d'sin(On’'(oi d*X H’2) 

It depends on two covariance matrices. The noise covariance matrix is known. It is 
a diagonal matrix with diagonal elements equal to q. The state covariance matrix 
on the other hand has to be computed. This was done through simulations of the 
system. Equation (4.2) can be regrouped into one matrix equation. 

£(r„'(OI = M„ W, 

As we want to optimize our relation for a zero failure of actuator 1 we have 
Ujg*(0 = 0 and our failure covariance matrix will be 

Mf = ^[u/(0Ui*l01 

The covariance ratio becomes 

(4.3) 

c = H- 

W J M„ W, 

ana 

Without sensor noise is simply 
M„ = C* £|X(0 X’’(0] C* 

With W C =0 the denominator is zero and c^ is infinite. Consequently in that 
case parity relations are the optimized relations. 
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Generally, this is not true and we have to find the vector that maximizes 
the ratio 


W, Mf 


( 4 . 4 ) 


W ' M„ W, 

ana 


Seen this way this doesn’t seem very easy. The idea is to first find a transformation 
for which the image of is the identity matrix. Let the eigenvalue decomposition 


of M„ be 
n 


Mjj = P D P’ with D a diagonal matrix: D = 


As the eigenvalues of D are positive, it can be written as 


D = D ^/2 1 ) 1/2 D ^/2 = 


(X,)l/2 




P is the transformation changing M into I. 


M = P D^/2 I (P 


Using this transformation the ratio can be rewritten as 

y^Tp jji/2 j)-l/2 pT w p D-1/2 1)1/2 pTyy 

a I a 

^TpDl/2 J D^/2p^W 

a a 

If we pose P^ and P^Mj P the ratio becomes 
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K Mft 




Finding the maximizing .this ratio is quite easy. If we impose that ||W^J1=1, 
then the denominator is one. The maximizing the numerator will be given by 
the eigenvector of correspoding to the largest eigenvalue. This technique is 
similar to what was done in 3.4. Once the best is found, we get the 

corresponding with 


W, = P d-'/2 w^, 

Optimized relations that will be generated here have the same number of 
coefficients, 147, as previous parity relations. They also have the norm of the 
sensor coefficients equal to one. This will enable us to compare the performances of 
the two types of relations. 


One single actuator 1 relation was generated for each of the two levels of 
noise. They were then tried in the two failures simulation. For the 10'^ level the 
relation was first tried without simulating any noise to see how the state affects the 
residual. The result is presented in figure 4-1. This particular case emphasizes the 
fact that letting through a little amount of the state can dramatically improve the 
failure signature. Figure 4-2 shows the simulation taking into account noise at the 
10'“* level. The equivalent simulations for the 10'^ level are shown without noise in 
figure 4-3 and with noise in figure 4-4. In this case much more of the state is visible 
in the residual. In return we should have a larger failure signature, as this relation 
must have a better covariance ratio than the preceding one for a level of noise of 
10'®. This cannot be seen in figure 4-3 which means that the increase must be quite 


small. 



Sensor 1 failure Actuator 1 failure 


Figure 4-2: Noisy residual of the model based optimized relation 
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Figure 4-4: Noisy residual of the 2“*^ model based optimized relation 
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4.3 Estimation Using the Optimization 

We now want to use the covariance ratio with the eigenvalue decomposition 
method. From section 3.4 we already know how to find the covariance in the no-fail 
case. The relation is of the form 

£[r2(0! = p;m„ P, 

We also know how to compute the covariance in the fail case. If the failure 
signature is W^Uj*(») 

e(r,2(0] = £:[r^2(0| + wl £:iu,'(0 u,*''(01 

And if our estimation uses n information vectors, we have 

Eiu,*(0 u/'to! =-•£ lUi'li) Ui*’(i)l 

" i=i 

The covariance ratio is 

Wj£(U,*(0Ui”(O| 
c, = 1 + 

py p, 

To be able to use the transformation technique of the previous section we want to 
express only as a function of P^.. As = f W^' , W^ ], if we define Mj. as 



Cj can be rewritten as 


This equation is equivalent to equation (4.3) except that now we are estimating all 
the coefficients of the relation instead of only the actuator ones. 

The optimized estimation was performed using 200 parity information 
vectors. First noise of 10"* was introduced. The resulting relation is shown in 
figure 4-5 for a simulation taking into account sensor noise. We somehow expected 
to get the same residual as in figure 4-2. Instead we have a totally different relation 
with a larger failure signature and a much higher level of residuals in the no fail 
case. The covariance ratio given by the estimation is 0.79 10^. Obviously, this 

A 

number does not. correspond to what we really have in the simulation. For the 10 
level of noise, the estimated relation is shown in figure 4-6. Here we have an even 
greater difference between the estimated ratio, 0.57 10^, and the simulation. This 
bad estimation of the covariance ratio indicates that the covariance matrices 
needed for the estimation cannot be reliably computed with 200 information 
vectors. Sensor noise is a likely source for the bias in the estimation. In section 3.4 
the same amount of information gave good results even when sensor noise was 
introduced. But if we compare figures 3-13 and 4-5, it seems that a small amount 
of information and sensor noise affects the estimation of optimized relations more 
than that of parity relations. 

In a second attempt the estimation was carried over 1000 time steps. Figure 
4-7 shows the result for the 10"* level of noise and figure 4-8 for the 10*^ level. The 
corresponding estimated covariance ratios are 0.18 10 and 0.13 10 . These 
numbers correspond more to what we actually see in the simulations. The most 
interesting point is that these relations have better performances than the ones of 



4-6 
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section 4.2. Previously the optimization was done only with respect to the actuator 
coefficients. Once these were chosen the rest of the relation was given by the 
matrix D*. This means that only 21 independent relations were considered. Now 
the optimization is done with respect to all the coefficients and we are comparing 
147 independent relations of which 21 are the previous ones. It turns out that the 
latter are not the best ones. 

Optimized relations are completely different from parity relations as they 
have a high level of residual in the no-fail case. As a consequence the effect of the 
sensor 1 failure is no longer visible in the residuals. This might also be the case for 
any actuator 1 failure other than the ’’zero” failure. 


ESTIMPTED R7 
. 00 
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Chapter 5 
Conclusion 

5.1 Summary 

In chapter 2, parity relations were generated using a model of a grid 
structure. The goal was to investigate the sensitivity of these relations to modeling 
errors. 5% error in the model parameters was enough to prevent us from detecting 
any actuator failure, and with 10%, no failure could be detected. As no accurate 
model can be found before the structure is built, chapter 3 focused on the 
estimation of parity relations using the actuator inputs and sensor outputs. 

Three on-line estimators were studied. First a Kalman filter was tried. The 
resulting estimated relation had a reduced detection efficiency compared to the 
relation generated using the model, due mainly to a reduced failure signature. This 
trend was observed in all estimators. Also this technique requires an amount of 
memory proportional to the square of the number of coefficients in the relation. As 
parity relations for large space structures will have many coefficients and memory 
on a space based computer is limited this can be a major problem. Consequently, a 
second estimator, for which the storage requirement is only proportional to the 
number of coefficients, was developed. It was based on the minimization of a 
structural distance. The results were comparable to what we had with the Kalman 
filter. However, the required estimation time is much longer. But this is not a 
problem in our case. This technique is the most interesting for single sensor parity 
relations. For a single actuator relation, however, a more accurate estimator is 
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required. The third technique used an eigenvalue decomposition of a covariance 
matrix. It proved to be slightly more accurate than the others with the same 
memory requirement as the Kalman filter. 

With a few modifications, the last technique was used to generate linear 
relations other than parity relations. Unlike the latter, these relations take into 
account the effect of the failure on the input-output relations. As a result they are 
optimized for the detection of a particular failure. The main interest of such 
relations is their ability to detect failures in very noisy conditions. It is likely that, 
in these conditions, they might only detect the failure they are designed for. 
However, components can usually fail in only one or a few different ways, so we can 
generally run a different relation for each different failure mode of the same 
component. If this b computationally too heavy, another solution is to modify the 
covariance ratio so as to optimize a relation for two 6r more failures.' For example, 
to detect the failure signatures Fl and F2, the covariance ratio would be 

1/2 covariance of Fl + 1/2 covariance of F2 

c^= 

covariance of the residual without failure 

5.2 Recommendations 

The last chapter showed the importance of the failure signature in the design 
of failure detection relations. Consequently, further research using this concept 
could be performed. In particular, the following points could be investigated. 

As a remedy to the complexity of the double eigenvalue decomposition, it 
would be interesting to use the covariance ratio concept with a simplified 
estimator. For example, we could define a structural distance as the inverse of the 
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covariance ratio and try to decrease it. But as this gives a nonlinear expression, the 
general framework of section 3.3 is not applicable and a new approach must be 
found. 

We have seen that estimated optimized relations have better performance 
than model based ones. But, whenever possible, it is more convienient to build 
relations with a model of the system rather than estimate them. Thus, further work 
could be done on the improvement of model based optimized relations. 

Finally, we can use the failure signature idea with parity relations. If we take 
the number S of time steps to be greater than the minimum, N-bl, we can choose 
between different independent relations. We can then choose the one with the best 
covariance ratio. For a parity relation, as the state is no longer present, this ratio is 

covariance of the failure signature 
^ ^ 

covariance of the residual due to noise 
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phi. GAMA. C AND FEEDBACK MATRICES OF THE DISCRETIZED SYSTEM 
PHI MATRIX 


0.97304D+00 0.98868D-01 0.000000+00 0.000000+00 O.OOOOCC-CO 
O.OOOOOD+00 O.OOOOOD+00 O.OOOOOD+00 0.000000+00 O.OOOOOD+00 
O.OOOOOD+00 0. 000000+00 0.000000+00 0.000000+00 0.000000+00 
0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 
-.536360+00 0.973090+00 0.000000+00 0.000000+00 0.000000+00 
0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 
0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 
0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 
0.000000+00 0.000000+00 0.919260+00 0.968940-01 0.000000+00 
0.000000+00 0.000000+00 O.OOOOOD+00 0.000000+00 0. 000000+00 
0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 
0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 
0.000000+00 0.000000+00 -.159040+01 0.919520+00 0.000000+00 
0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 
0.000000+00 0.000000+00 O.OOOOOD+OO 0.000000+00 0.000000+00 
0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 
0.000000+00 0.000000+00 0.000000+00 0.000000+00 0.620950+00 
0.861530-01 0.000000+00 0.000000+00 0.000000+00 0.000000+00 
0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 
0.000000+00 0.000000+00 0.000000+00 O.OOOOOD+OO 0.000000+00 
0.000000+00 0.000000+00 0.000000+00 0.000000+00 -.703610+01 
0.623450+00 0.000000+00 0.000000+00 0.000000+00 O.OOOOOD+OO 
0.000000+00 O.OOOOOD+OO 0.000000+00 O.OOOOOD+OO O.OOOOOD+OO 
O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO 
0.000000+00 O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO 
O.OOOOOD+OO 0.183130+00 0.696390-01 O.OOOOOD+OO O.OOOOOD+OO 
O.OOOOOD+OO 0.000000+00 O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO 
O.OOOOOD+00 O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO 
O.OOOOOD+OO 0.000000+00 O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO 
O.OOOOOD+OO -.135340+02 0.191590+00 O.OOOOOD+OO O.OOOOOD+OO 
0.000000+00 O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO 
O.OOOOOD+Ob 0.000000+00 0.000000+00 O.OOOOOD+OO 0.000000+00 
O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO 
0.000000+00 O.OOOOOD+OO O.OOOOOD+OO -.353180+00 0.467540-01 
0.000000+00 O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO 
O.OOOOOD+OO 0.000000+00 O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO 
O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO 0.000000+00 O.OOOOOD+OO 
O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO -.177650+02 -.332420+00 
O.OOOOOD+OO 0.000000+00 O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO 
O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO 
0.000000+00 O.OOOOOD+OO 0.000000+00 O.OOOOOD+OO O.OOOOOD+OO 
O.OOOOOD+OO 0.000000+00 O.OOOOOD+OO 0.000000+00 O.OOOOOD+OO 
-.958840+00 O. 469220-02 O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO 
0.000000+00 0.000000+00 O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO 
0.000000+00 0.000000+00 0.000000+00 0.000000+00 O.OOOOOD+OO 
0.000000+00 O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO 
-.421280+01 -.901720+00 O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO 
O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO 
0.000000+00 O.OOOOOD+OO 0.000000+00 O.OOOOOD+OO O.OOOOOD+OO 
0.000000+00 O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO O.OOOOOD+OO 
0.000000+00 O.OOOOOD+OO -.867910+00 -.119250-01 0.000000+00 
0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 
0.000000+00 0.000000+00 0.000000+00 0.000000+00 O.OOOOOD+OO 
0.000000+00 0.000000+00 O.OOOOOD+OO 0.000000+00 O.OOOOOD+OO 
O.OOOOOD+OO 0.000000+00 0.15485D+02 -.787240+00 0.000000+00 
0.000000+00 0.000000+00 O.OOOOOD+OO 0.000000+00 O.OOOOOD+OO 
O.OOOOOD+OO 0.000000+00 0.000000+00 0.000000+00 O.OOOOOD+OO 
0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 
O.OOOOOD+OO 0.000000+00 O.OOOOOD+OO 0.000000+00 -.810230+00 
-.142190-01 0.000000+00 O.OOOOOD+OO O.OOOOOD+OO 0.000000+00 
0.000000+00 O.OOOOOD+OO 0.000000+00 O.OOOOOD+OO O.OOOOOD+OO 
0.000000+00 0.000000+00 O.OOOOOD+OO 0.000000+00 O.OOOOOD+OO 
0.000000+00 0.000000+00 0.000000+00 0.000000+00 0.197150+02 
-.725170+00 0.000000+00 O.OOOOOD+OO 0.000000+00 0.000000+00 
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0 . 000000+00 O. 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 
0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 
0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 
0.000000+00 -.507000-01 -.204210-01 0.000000+00 0.000000+00 
0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 
0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 
0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 
0.000000+00 0.445250+02 0.617610-01 0.000000+00 0.000000+00 
0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 
0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 
0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 
0.000000+00 0.000000+00 0.000000+00 0.915790+00 0.319060-02 
0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 
0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 
0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 0 . 000000+00 
0.000000+00 0.000000+00 0.000000+00 -.135090+02 0.104180+01 


GAMA MATRIX 


- .746900-06 
- . 148590-04 
0.292700-03 
0.576580-02 
- .468490-05 
- .869620-04 
0.209250-05 
0.346700-04 
- .272180-03 
-.357340-02 
- . 169030-04 
-.363540-04 
0.494850-05 
-.410230-04 
- . 169310-03 
O. 184400-02 
- .458060-04 
O. 194110-02 
0.417320-05 
0.669480-03 


O. 104710-03 
0.208310-02 
O. 174380-03 
0.343510-02 
-.235670-03 
-.437460-02 
-.396230-03 
-.656490-02 
O. 162250-03 
0.213010-02 
0.203690-03 
0.438070-03 
-.750410-04 
0.622090-03 
0.621370-04 
- .676720-03 
0.738970-04 
- .313150-02 
- .'238470-05 
- .382550-03 


- . 105590-03 
-.210060-02 
O. 179580-03 
O. 353750-02 
0.229320-03 
0.425670-02 
0.397900-03 
0.659250-02 
O. 168010-03 
0.220570-02 
-. 186770-03 
-.401680-03 
0.718160-04 
-.595360-03 
0.883690-04 
- .962420-03 
0.724320-04 
-.306940-02 
- . 226770-05 
-.363790-03 


-. 147590-03 
-.293610-02 
0.362520-04 
0.714120-03 
0.299550-03 
0.556040-02 
0.365790-04 
0.606060-03 
-.222380-03 
-.291950-02 
O. 137530-03 
0.295780-03 
0.631410-04 
-.523440-03 
-. 189380-04 
0.206250-03 
0.251560-04 
- . 106600-02 
0.248890-05 
0.399270-03 


O. 147430-03 
0.293300-02 
0.292210-04 
0.575620-03 
-.301340-03 
-.559350-02 
- . 347190-04 
-.575230-03 
-.218570-03 
-.286950-02 
-. 142240-03 
-.305900-03 
- .621950-04 
0.515600-03 
-.216040-05 
0.235290-04 
0.234760-04 
- .994830-03 
0.254460-05 
0.408200-03 


- .616150-06 
- . 122580-04 
0.233390-03 
0.459760-02 
-.252290-05 
- .468320-04 
-.701750-06 
- . 116270-04 
O. 191890-03 
0.251920-02 
-.214900-04 
-.462180-04 
0.243070-05 
-.201500-04 
-.210260-03 
0.228990-02 
-. 132640-04 
0.562100-03 
-.321420-05 
-.515630-03 


C MATRIX 


-.297440+01 0.000000+00 0.363520+01 0.000000+00 0.264480+01 
0.000000+00 -.306310+01 0.000000+00 -.439700+01 0.000000+00 
-.251550+01 0.000000+00 -.532520+01 0.000000+00 0.448740+01 
0.000000+00 -.351850+01 0.000000+00 -.551090+01 0.000000+00 
-.153290+01 0.000000+000.319720+01 0.000000+00 -.109100+01 
0.000000+00 -.374940+01 0.000000+00 0.146060+01 0.000000+00 
0.192900+01 0.000000+00 0.976480+00 0.000000+00 0.242190+01 
0.000000+00 0.381260+01 0.000000+00 0.572960+00 0.000000+00 
-.457930+00 0.000000+00 -.190700-01 0.000000+00 -.349540+01 
0.000000+00 0.168910+01 0.000000+00 0.717730-01 0.000000+00 
-.414990+01 0.000000+00 -.334710+01 0.000000+00 0.182330+00 
0.000000+00 -.155360-01 0.000000+00 -.119650+00 0.000000+00 
0.000000+00 0.210700-01 0.000000+00 0.354520-01 0.000000+00 
-.507770-01 0.000000+00 -.942700-01 0.000000+00 0.455600-01 
0.000000+00 0.933610-01 0.000000+00 -.521660-01 0.000000+00 
0.475920-01 0.000000+00 0.153350+00 0.000000+00 -.119900+00 
0.000000+00 O. 296660-01 0.000000+00 0.594070-02 0.000000+00 
-.649260-01 0.000000+00 -.826020-02 0.000000+00 -.613740-01 
0. 000000+00 -.651930-01 0.000000+00 -.432360-01 0.000000+00 
-.165470-02 0.000000+00 0.487170-01 O. 000000+00 0.127940+00 
0.000000+00 -.123980-03 0.000000+00 0.474500-01 0.000000+00 
-.543590-03 0.000000+00 -.166960-03 0.000000+00 0.538810-01 
0.000000+00 -.984980-02 0.000000+00 0.168970-02 0.000000+00 
-.161040+00 0. 000000+00 -.275260-01 0.000000+00 -.161610+00 
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feedback MATRIX 


O. 14950D+03 O. 13786D+02 -.726200+02 -.222250+02 -.205480+03 
-.223940+02 -.178270+03 0.377390+01 0.190740+03 0.282130+02 
-.258060+03 0.157410+02 -.369490+03 -.994980+01 0.122440+04 
-.253870+02 -.137040+04 -.237780+02 -.546010+03 0.115800+03 
0.502900+02 0.435230+01 -.849090+02 -.114590+02 -.146380+02 
0.573190+00 0.244510+02 0.784550+01 0.139880+02 -.193590+00 
0.349320+02 -.283410+01 0.196250+02 -.612470+000.457360+03 
-.958770+01 0.770380+02 0.292860+01 -.636170+02 0.755710+01 
0.116230+03 0.150960+02 -.104810+03 -.134650+02 0.598820+02 
-.267820+01 -.628540+01 0.230820+01 0.131510+03 0.119160+02 
-.212460+03 0.964740+01 -.357540+03 -.607710+00 0.519460+03 
-.471520+01 0.355460+03 -.612230+01 -.286550+03 0.547490+02 
0.146200+03 0.233740+02 -.312850+03 -.612440+01 0.338970+03 
0.165420+02 0.232680+03 0.163690+02 -.206040+03 -.147250+02 
0.266660+02 0.422900+01 -.237550+03 0.109100+02 -.197950+04 
0.919900+02 0.160740+04 0.200360+02 0.129760+04 -.116550+03 
0.143800+03 0.215570+02 -.292160+03 -.472200+01 0.428630+03 
0.342920+02 0.222860+03 0.105100+02 -.242840+03 -.184680+02 
-.126860+03 0.532270+00 -.276470+03 0.836690+01 -.190630+04 
0.870110+02 0.152430+04 0.222540+02 0.136260+04 -.129810+03 
0.261480+03 0.413270+02 -.427510+03 -.136640+02 0.599540+03 
0.321570+02 0.105170+03 0.157660+02 0.528980+02 -.679320+01 
-.211830+03 0.346040+02 -.666170+03 0.110610+02 -.289970+04 
0.138200+03 0.869760+03 0.184810+02 0.720960+03 -.558490+02 


